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CA 

Clear  aperture 

ESI' 

Electrostatic  Unit 

FWHM 

Full  width  at  half  maximum  Intensity 

GDL 

Glass  development  laser 

Np 

Neper 

ns 
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PS 
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Synthesized  Neodymium  Amplifier  Analysis  Program 
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A Preliminary  Analysis  of  GDI.  Performance 
for  Short  and  I-ong  Pulse  Operation 


The  development  of  a three-dimensional  gain  c«Kle  for  system  performance  predictions 
has  been  accomplished.  This  code.  Synthesized  Neodymium  Amplifier  Analysis  Program 
(SNAAP),  has  been  used  to  predict  the  operating  performance  of  the  University  of  Rochester's 
Glass  Development  Laser  (GDL).  In  its  present  form,  the  code  essentially  traces  the  propa- 
gation of  a laser  pulse  represented  by  a three-dimensional  matrix  through  a given  optical 
system. 

For  the  GDI.  system  a pulse  whose  spatial  and  temporal  waveforms  were  initially 
gaussian  was  convolved  with  a spatially  apodizing  function  and  then  propagated  through  the 
system.  The  pulse  is  represented  by  a matrix  (m,  n,  j).  Each  m1*1  x n1*1  x j1*1  beam  element 
represents  an  intensity  (w/cm2).  The  code  then  propagates  each  such  element  and  o|*?rates 
on  it  according  to  the  type  and  sequence  of  components.  Physically  there  are  three  types  of 
optical  components: 

1.  Amplifiers  which  include  all  gain  media  where  passive  and  reflective 
losses  are  calculated  where  indicated. 

2.  Spatial  filters  where  an  assumed  transfer  function  operates  on  each 
beam  element  and  readjusts  its  intensity  accordingly. 

3.  Attenuators  which  include  all  other  optical  elements  such  as  lenses, 
polarizers,  Pockels'  cells,  etc. 

In  its  present  form,  most  calculated  parameters  including  B-integral  (both  total  and 
delta  B per  stage),  each  intensity  (m1*1  x n1*1  x j1*1),  etc.  are  stored  in  arrays  for  future 
calculations  or  manipulation.  A great  deal  of  flexibility  and  universality  has  been  built 
into  the  code  to  enable  both  present  and  future  users  ease  of  modification  and  an  under- 
standing of  the  codes  operation. 

The  input  parameters  used  for  this  preliminary  analysis  appear  in  Table  1.  The 
stored  energy  density  profiles  for  each  amplifier  are  input  in  the  format 

Eg(r')  = + E2  r,P  (1) 

The  specific  gain  for  several  types  of  Nd:doped  glass  (LHG-7,  Q-88,  ED-2,  EV-2)  are 
implemented  in  the  code  and  the  appropriate  values  used  depending  on  the  glass  type  speci- 
fied as  input. 
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Presently  the  spatial  filter  transmission  function  used  for  the  pinhole  calculation  is 
taken  from  the  Lawrence  Livermore  laboratory  1975  Annual  Report.  The  function  used  is 

T(B>  - exp  [-(B/3.4)8]  (2) 

Associated  with  each  mth  x 0th  x jth  intensity  element  Is  a B-integral  value  B(m,  n,  j)  which 
attenuates  each  intensity  element  at  a pinhole  according  to  Equation  (2).  Modifications  to 
Equation  (2)  are  presently  being  discussed  and  may  be  implemented  into  the  code  at  a later 
date. 

The  effect  of  circular  polarization  is  not  at  present  contained  within  the  code  and  is 
handled  by  inputting  n.,/1.5  in  place  of  the  n 2 values  listed  in  Table  1.  Another  parameter 
which  is  not  presently  implemented  into  the  code  is  fill  factor.  For  the  purpose  of  sim- 
plicity, the  code  exactly  fills  the  aperture  of  an  optical  element  with  mxn  spatial  elements 
of  specified  profile.  As  an  example,  suppose  we  have  a 15-mm  diameter  clear  aperture 
amplifier,  area  element  mxn  (actually  intensity)  is  magnified  by  a factor  of  four  but  occupies 
the  same  relative  position  in  amplifier  II  as  it  did  in  amplifier  I (center  transforms  to 
center,  and  edge-to-edge,  etc.). 

In  the  beginning  of  this  report  it  was  mentioned  that  the  spatial  profile  of  the  beam  was 
a gaussian  convoluted  with  a supergaussian.  More  specifically,  the  actual  profile  is 

O IQ 

I(r)  = exp(-(r/0.3)  -(r/0.53)  ] (3) 

where  r is  in  cm.  This  spatial  profile  is  the  input  to  the  first  (16-mm)  amplifier.  Since 
the  fill  factor  is  0.  26  and  the  code  does  not  automatically  adjust  for  vacuum  spatial  filter 
magnification  (it  adjusts  for  differences  in  clear  aperture  (CA>  only),  it  is  necessary  to 
correct  for  the  fact  that  the  first  two  amplifiers  (i.e.  16-mm  and  30-mm)  are  underfilled. 
Basically,  we  wish  to  map  the  0.001  I beam  edge  onto  the  edge  of  the  40-mm  (third)  ampli- 
fier. From  the  magnification  factors  listed  in  Table  1 this  is  accomplished  by  a simple 
linear  transformation 

16  mm—  12.60  mm  : rj,  = 6.3  mm  ^ 

30  mm  — 28.99  mm  : r^  14.5  mm 

We  designate  these  new  diameters  as  the  effective  clear  aperture  for  the  respective  ampli- 
fiers. The  radial  gain  profiles  for  the  code  are  all  normalized  to  r/r  , where  r goes  from 
0 to  rQ.  In  order  to  correctly  correlate  the  initial  beam  profile  with  the  amplifier  gain 
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profile,  It  is  necessary  to  transform  the  initial  stored  energy  density  profiles  such  that 
from  Equation  (1)  we  have 

E’  E2(0.79)1,75  = 0.16  (5|* 

for  the  16-mm  amplifier,  and 

E!,  Eo<0.97>2,87  0.12  (6 1* 

for  the  30-mm  amplifier. 

Besides  the  system  output  parameters  provided  by  the  code,  numerous  component 
operating  parameters  are  provided  at  each  optical  interface.  Due  to  the  large  number  of 
system  components,  the  latter  are  provided  mainly  as  a diagnostic  feature.  Once  the  system 
components  are  specified,  the  basic  input  to  the  code  is  the  initial  beam  intensity.  The  code 
calculates  numerically. 

Eo  = ///  V*"’  l)  dr  d0  dt 

where 

Eq  = Output  pulse  energy 

1.  Input  intensity  profile 

The  first  code  runs  for  the  GDL  system  were  at  a temporal  pulse  length  of  50  ps 
(FVVHM).  The  test  cases  run  were  for  inputs  to  the  first  (16-mm  diameter)  amplifier  of 
5 pJ/cm2  (1.3  pJ  total  beam  energy)  and  6 pJ/cm2  (1.6  pJ).  Figures  1 and  2 are  CalComp 
plots  of  the  output  intensity  profiles  after  SF-7  (see  Table  1)  for  these  two  cases,  res)>ec- 
tively.  These  graphs  are  the  result  of  taking  a series  of  21  radial  increments  across  the 
beam,  propagating  these  through  the  system  and  adjusting  the  Intensity  of  each  m1*1  x n1*1  x 
jth  element  for  small-scale  B-integral  losses  through  the  various  vacuum  spatial  filters  in 
the  system.  Figure  2 is  a rather  dramatic  representation  of  the  effects  of  B-integral  on  a 
beam  focused  through  the  pinhole  of  a vacuum  spatial  filter.  Figure  3 shows  a plot  of  the 

o 

intensity  profile  for  the  same  6 pJ/cm  input  as  Figure  2 but  taken  after  the  90-mm  amplifier. 
This  is  a representation  of  the  near  field  of  the  beam.  One  can  also  see  the  spatial  and 
temporal  evolution  of  the  beam  (for  the  6 pJ/cm2  input  case)  from  Figures  4 through  8 
inclusive).  These  are  plots  of  the  intensity  profiles  for  the  40-mm,  first  64-mm,  second 


♦The  primes  denote  the  effective  values  of  E„  from  E (r/r  ) E,  + 
Equation  (1).  * 8 ° 


E'2(r/r0)> 


see 
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64-mm,  third  64-mm  and  the  90-mm  amplifier,  respectively.  The  radial  stored  energy 
density  profiles  for  these  amplifiers  account  for  the  higher  intensities  on  the  edge  of  the 
beam,  f igures  9 and  10  are  plots  of  the  output  energy  densities  as  a function  of  normalized 
radial  beam  position  for  the  5 pj/cm2  and  6 pJ/cm2  runs,  respectively.  Again  the  outputs 
plotted  are  those  through  SF-7  (see  Table  1).  As  mentioned  before,  associated  with  each 
beam  intensity  element  is  a B-integral  element  B(m,  n,  j».  Figures  11  and  12  are  (dots  of 
the  system  B-integral  values  for  these  two  runs,  respectively. 

The  final  GDL  system  output  energy  is  obtained  by  integrating  the  beam  intensity  pro- 
file in  both  space  and  time.  This  is  accomplished  by  integrating  the  functions  appearing  in 
Figures  9 and  10.  After  performing  the  required  integration,  values  of  40.  1 J for  the 
5 pJ/cm2  case  and  43. 1 J for  the  6 pj/cm2  case  were  obtained.  Due  to  the  long  computing 
time  necessary  to  run  test  cases  of  21  spatial  points  across  the  beam,  it  was  decided  that 
for  sufficient  accuracies  four  points  across  the  beam  would  be  adequate  to  predict  system 
performance.  These  points  are  0,  0.35,  0.8,  and  1.0  r/rQ.  Note  from  Figure  9 that  the 
l>eak  occurs  at  the  0.  8 r/r()  point.  The  difference  in  integrations  of  4 (joints  vs  21  points  is 
4.7*2.  For  the  6 pJ/cm2  case,  the  difference  is  8.9*2;  however,  in  this  case  the  maximum 
does  not  occur  at  the  0.  8 r/r  radial  point  due  to  B-integral  considerations.  In  any  event, 
the  4-point  integration  produces  more  pessimistic  results  than  the  21-point  case. 

In  Table  2,  the  output  energy  vs  input  energy  for  the  four  temporal  cases  run  (viz. 

50  ps,  100  ps,  200  ps . and  1 ns)  is  presented. 

Figures  13  through  21  inclusive  present  the  beam  output  energy  density  through  SF-7 
(see  Table  1)  as  a function  of  radius  for  the  nine  input  energies  at  50  ps  listed  in  Table  2. 

The  vertical  axis  is  in  J/cm  and  the  beam  radius  is  normalized  to  one.  Table  2 presents 
the  input  in  J whereas  in  several  places,  such  as  the  plots  mentioned  above,  the  input  is  given 
in  J/cm2.  The  assumed  conversion  for  the  two  has  been  calculated  to  be  0.  2626  cm"  (the 
beam  effective  aperture  area  at  the  input  to  the  16-mm  diameter  amplifier).  The  plots  in 
f'igures  13  through  21  are  in  descending  order  of  output  energy  (or  input  energy  density). 

Figures  13  through  15  clearly  show  the  effects  of  B-integral  through  SF-7  (see  Table  1). 

2 2 

The  peak  output  energy  density  for  the  6. 5 pJ/cm  input  case  is  0.  82  J/ cm"  at  the  edge 

2 2 
(0.  8 r/rQ);  however,  for  the  5 pJ/cm  case,  the  output  energy  at  0.  8 r/rQ  is  0.  88  J/cm  . 

Figure  22  is  a plot  of  the  output  energy  as  a function  of  input  energy  for  the  nine  cases  run 

at  50  ps.  Similar  data  for  the  other  three  temporal  runs  is  presented  (see  Table  2).  Only 

brief  mention  of  some  interesting  points  of  these  runs  will  be  made  since  many  of  the  above 

comments  apply  for  all  four  temporal  pulse  lengths.  Because  of  the  effects  of  the  vacuum 
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spatiiil  filters,  the  output  radial  beam  profile  (at  100  ps)  is  almost  flat  for  the  20  pj/cm2 
input  case  (see  Figure  23).  This  is  due  to  the  higher  intensity  at  the  0.8  r/r  ( beam  radius. 
Figure  33  is  a plot  of  the  temporal  evolution  of  this  pulse  at  the  0.  8 r/r()  |>oint.  The  severe 
dip  in  the  center  of  the  pulse  occurs  due  to  focusing  through  SF-7.  Figure  32  is  a plot  of  the 
output  energy  vs  input  energy  for  the  100-ps  case  (see  Table  2).  Figures  34  through  44 
inclusive  contain  plots  of  the  output  energy  density  profiles  of  the  11  cases  at  200  ps  listed 
in  Table  2,  with  Figure  45  showing  the  output  vs  input  energy  curve  for  these  cases.  Finally, 
11  cases  for  an  input  pulse  width  of  1 ns  were  run.  The  output  energy  density  (through  SF-7) 
for  these  cases  is  plotted  in  Figures  46  through  56  inclusive.  For  this  case,  the  output  is 
saturation  flux  limited  and  not  B-integral  limited  as  in  the  previous  three  cases.  As  an 
example,  the  highest  input  energy,  1000  pj/cm2  (263  pJ),  has  a AB  through  the  last 
amplifier  stage  (SF-7)  of  1.37.  However,  the  assumed  flux  limit  of  4 J/cnr  is  exceeded  at 
the  output  of  the  second  64-mm  rod  amplifier.  Assuming  a system  could  tolerate  a higher 
flux  than  4 J/cm2  at  any  surface  (say,  for  example,  6 J/cm2),  it  is  interesting  to  note  that 
the  gain  of  the  second  64-mm  amplifier  is  down  by  a factor  of  four  from  the  small  signal 
gain.  At  such  high  fluxes,  the  concept  of  gain  reduces  to  cot  sidering  an  amplifier  as  an 
additive  medium,  rather  than  a gain  medium;  i.e., 

high  flux 

E = GE  - E = E + A E 

° limit 

Although  both  of  the  above  expressions  are  appropriate  in  either  case,  the  latter  is  more 
useful  at  high  incident  flux  levels.  In  order  to  insure  that  no  system  component  exceeds  the 
flux  limit  (4  J/cm2)  the  input  energy  density  must  not  exceed  145  pJ/cm2.  At  an  input  of 
100  pj/cm2  (26  pJ)  the  90-mm  amplifier's  gain  is  down  by  almost  a factor  of  three  from  the 
small  signal  gain.  This  is  pointed  out  to  illustrate  the  fact  that  the  GDL  system  as  presently 
configured  is  saturatiop  flux  limited  under  long  pulse  operation  1 ns)  at  low  input  energies 
(^1  pJ  input).  This  point  will  not  be  dwelt  upon  here  but  the  results  of  saturation  will  be 
presented  graphically  later  on  in  this  report. 

For  each  of  the  four  temporal  pulse  length  cases  analyzed  (50  ps,  100  ps,  200  ps,  and 
1 ns)  input  vs  output  energy  is  plotted  (Figures  58  through  61).  Each  graph  also  shows  the 
total  B-integral  values  (Np)  for  the  system  (at  the  output  of  SF-7),  as  well  as  the  delta 
B-integral  through  the  last  stage  (SF-6  to  SF-7  "pinhole  to  pinhole";  see  Table  1).  In  each 
case  values  of  B-integral  vs  input  energy  for  the  center  (0.0  r/rQ)  of  the  pulse  and  the  edge 
(0.  8 r/rQ)  have  been  plotted.  The  latter  in  all  cases  is  the  higher  of  the  two  plots.  Finally, 
the  plots  have  been  taken  and  system  output  energy  vs  AB  (last  stage)  = 2.5  has  been  cross 
correlated.  This  value  was  chosen  conservatively  and  is  an  indication  of  system 
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performance  at  this  given  B-integral  limit.  It  is  interesting  to  note  (although  it  is  not 
obvious  from  the  graph)  that  for  the  200-ps  case  (Figure  60),  the  system  total  B-integral 
CuB  on  the  plot  (center  vs  edge))  crosses  over  for  high  values  of  input  energy,  typically  at 
L'jn  • 16  pJ  where  the  2JB  is  the  same  for  both  center  and  edge.  This  crossover  is  not  shown 
in  Figure  60  since  the  plot  does  not  extend  beyond  a 12  p.J  input  energy.  Figure  61  is  again 
similar  to  the  previous  three;  however,  in  this  case  the  flux  (J/cm'S  rather  than  the 
B-integral  for  both  the  center  and  the  edge  points  has  been  plotted.  This  is  the  limiting  case 
for  long  pulse  operation  of  GDL  since  the  B-integral  values  are  intensity-de pendent;  and  due 
to  saturation  at  long  temporal  pulse  length,  the  intensity  is  considerably  reduced  from  that 
of  the  shorter  pulses.  This  is  evident  in  Figure  62.  Input  energy  with  varying  scales  of  p.J 
corresponding  to  the  differences  in  pulse  length  is  plotted  vs  output  |»eak  power.  By  appro- 
priate choice  of  scale  the  vertical  axis  corresponds  to  constant  output  |>ower  at  any  of  the 
pulse  lengths  (full  scale  represents  0.8  TW).  The  dashed  line  on  this  graph  represents  the 
small  signal  gain  of  the  system.  Although  both  B-integral  and  saturation  effects  are  com- 
peting for  depleting  constant  system  output  gain  for  the  50-  , 100-  and  200-ps  cases,  it  is 
i still  obvious  that  saturation  effects  account  for  a |x>rtion  of  the  lO'/f  loss  in  output  power  for 

the  100-ps  case  over  the  50-ps  case  (at  the  peak)  and  the  approximate  20'J  loss  for  the  200-ps 
case  over  the  100-ps  case.  The  saturation  effects  for  the  1 ns  case  are  quite  obvious.  Here 
the  B-integral  effects  are  negligible  since  even  at  the  highest  input  energy  case,  1000  p.J  cm“ 
(263  pj).  the  total  system  B-integral  is  only  4. 1 with  a AB  (all  B-integrals  are  maximum 
spatial  and  temporal)  through  the  last  stage  (SF-7)  of  1.4. 

Figure  63  is  similar  to  Figure  62  except  here  AB  = 2.  5 has  been  plotted  through  the 
last  stage  (ABg())  for  the  three  shortest  temporal  pulse  length  cases  (50  ps,  100  ps,  and 
200  ps).  In  terms  of  the  1-ns  run  the  flux  limit  of  4 J/cm2  is  indicated.  Beyond  this  (joint, 
saturation  effects  are  so  limiting  that  even  if  new  substrates  withstanding  higher  flux  were 
found  very  little  additional  output  could  be  derived  from  this  system.  To  illustrate  this  more 
graphically  Figure  64  is  presented,  which  has  a compressed  horizontal  axis.  Here  it  is 
shown  for  the  aforementioned  run  that  at  an  input  of  50  p J/cm  a fourfold  increase  in  input 
energy  yields  an  approximate  15%  increase  in  output  power. 


One  final  run  for  the  GDL  system  will  be  made  at  500  ps  for  the  present  operating 
parameters  (Table  1).  This  run  should  supply  the  data  necessary  to  establish  the  following: 

1.  Is  the  AB90  = 2.5  cross  plot  of  system  performance  (Figure  63)  a 
straight  line  ? 

2.  Where  is  the  crossover  (joint  for  flux  limit  vs  B-integral  limit? 

When  this  run  is  made,  the  results  will  be  issued  as  an  addendum  to  this  report. 
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In  Figure  65  plots  of  the  output  of  GDI,  prior  to  SF-7  have  been  added  to  the  previous 
information  contained  in  Figure  63  in  order  to  illustrate  the  effects  of  B-integral  on  the 
system. 

The  overall  system  performance  is  given  in  Figure  66,  where  on-target*  output  energy 
and  output  peak  power  are  plotted  as  a function  of  input  pulse  length  full  width  half  maximum 
(FWHM).  The  output  energy  increases  linearly  below  200  ps,  with  a maximum  of  200  J at 
1 ns  limited  by  the  saturation  gain  of  the  amplifiers.  The  maximum  peak  power  of  the  system 
is  0.  7 TW  for  short  pulse  lengths  of  50  ps  or  less.  This  peak  power  is  limited  by  the  action 
of  the  final  spatial  filter. 

In  summary  the  GDI.  system  as  configured  in  Table  1 <one  arm  of  OMEGA-10)  will 
produce  the  following  outputs  (see  Figure  66): 

0.7  TVV  at  50  ps 
0.65  TW  at  100  ps 
0.6  TW  at  200  ps 
200  J at  1 ns 

Consequently,  OMEGA- 10  up  to  and  including  the  90-mm  diameter  amplifiers  will  produce: 

16.8  TW  at  50  ps 
15.6  TW  at  100  ps 
14.4  TW  at  200  ps 
4.  8 kJ  at  1 ns 


*Focusable  within  500  prad. 
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OUTPUT  INTENSITY  PROFILES 
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Fig.  1.  After  SF-7,  5 pJ/cm  (1.3pJ) 


Fig.  2.  After  SF-7,  6 pJ/cm  (1.6  pJ) 
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Fig.  3.  After  90-mm  Amplifier,  6 pJ/cni  Fig.  4.  After  40-mm  Amplifier,  6 pJ/cm“ 


MAX  P VALI1  * I lO1* 


X' 


MAX  n«im 

o <xa<jocn 


»' 


Fig.  5.  After  First  64-mm  Amplifier, 
6 pJ/cm2 


Fig.  6.  After  Second  64-mm  Amplifier, 
6 pJ/cm2 
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Fig.  7.  After  Third  64-mm  Amplifier, 
6 pJ/cm2 
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Fig.  8.  After  90-mm  Amplifier, 
6 pJ/cm2 
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OUTPUT  ENERGY  DENSITIES  AS  A FUNCTION  OF 
RADIUS  THROUGH  SF-7 
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Fig.  9.  5 pj/cm2 


Fig.  10.  6 pJ/cm2 


SYSTEM  B-INTEGRALS 


Fig.  11.  5 pj/cm2  Fig.  12.  6 pJ/cm2 


OUTPUT  ENERGY  DENSITIES  AS  A FUNCTION 
OF  RADIUS  AT  50  ps 
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Fig.  13.  6 . 5 pj/cm2  — 39 . 19  J Fig.  14.  6.0  pj/cm2~39.61  J 


Fig.  15.  5.  0 pj/cm2  — 38.  28  J 


9 


OUTPUT  ENERGY  DENSITIES  AS  A FUNCTION 
OF  RADIUS  AT  50  ps  (Cont) 
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Fig.  19.  2.0  pJ/cm2  — 19. 13  J 
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Fig.  20.  1.0  pJ/cm2  — 10.20  J 
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Fig.  21.  0.5  pJ/cm2  — 5.28  J 
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Fig.  22.  Output  vs  Input  Energy  at  50  ps 


OUTPUT  ENERGY  DENSITIES  AS  A FUNCTION 
OF  RADIUS  AT  100  ps 
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Fig.  23.  20.0  pj/cm2—  73.98  J 


Fig.  24.  18. OpJ/cm2—  74.75  J 
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Fig.  25.  16.0  pj/cm2  — 74  . 24  J 
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Fig.  27.  9.0  pj/cm2— 61. 12  J 


Fig.  26.  12.0  pJ/cm2—  70.90  J 
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Fig.  28.  7.0  p J/cm2  — 51.  87  J 
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Fig.  29.  6.0  pj/cm2  — 46  51 J 
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Fig.  30.  5.  0 jx J/ cm2 -*  40. 64  J Fig.  31.  4. 0 pJ/cm2- 34. 19  J 


ii 


Fig.  32.  Output  vs  Input  Energy  Fig.  33.  Temporal  Evolution  of  CXitput 

at  100  ps  Radial  Beam  at  0.8  r/r0 

OUTPUT  ENERGY  DENSITIES  AS  A FUNCTION 
OF  RADIUS  AT  200  ps 
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Fig.  34.  80.0  p«J/cm2  — 144.77  J 
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Fig.  35.  70. 0 pj/cm2  — 141.  84  J 
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Fig.  36.  60.0pJ/cm2  — 139. 12  J 
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Fig.  37.  52.  0 pj/cm2  — 136.47  J 
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Fig.  38.  48.0pJ/cm2  — 134.65  J 


Fig.  39.  44.  0 pj/cm2— 132.  28  J 
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OUTPUT  ENERGY  DENSITIES  AS  A FUNCTION 
OF  RADIUS  AT  200  ps  (Cont) 
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Fig.  40.  40.  0 n J/cm2  — 129.  20  J 


; 

4 


rr 

. 


f r TTTTl  "Pd 

:::::::  □ 


Fig.  41.  36.0pJ/cm2  — 125.  28  J 
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Fig.  42.  30.0  pj/cm2  — 117 
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Fig.  43.  20.  0pJ/cm2  — 98.01J 
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Fig.  44.  5.  0 p J/cm2  — 40.65  J 
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Fig.  45.  Output  vs  Input  at  200  ps 
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OUTPUT  ENERGY  DENSITIES  AS  A FUNCTION 
OF  RADIUS  AT  1 ns 


Fig.  46.  1000.0  n<J/cm2—  272.96  J Fig.  47.  850. 0 nJ/cm2  — 267.  77  J 
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Fig.  48.  650.0  nJ/cm2—  258.76  J 
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Fig.  49.  500.0  fiJ/cm2  — 249.43  J 


OUTPUT  ENERGY  DENSITIES  AS  A FUNCTION 
OF  RADIUS  AT  1 ns  (Cont) 
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Fig.  52.  200.0  pj/cm2  — 212. 66  J 
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Fig.  53.  100.0  pj/cm2  — 180.63  J 


Fig.  54  . 50.  0 pJ/cm2-»  145.62  J 


Fig.  55.  10.0  pj/cm2  — 66. 02  J Fig.  56.  5.0  pj/cm2  — 40.65 
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Fig.  57.  Output  vs  Input  Energy  at  1 ns 
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OUTPUT 
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Fig.  62.  Output  Power  vs  Input 
Energy 


Fig.  63.  Output  Power  vs  Input  Energy 
and  AB90  = 2.  5 
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Fig.  64.  Output  Power  vs  Input  Energy, 
Compressed  Scale 


Fig.  65.  Output  Power  vs  Input  Energy 
Showing  Effects  of  B-lntegral 
on  the  System 


Fig.  66.  System  Output  Performance 
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TABLE  1.  IN  PUT  PARAMETERS 


4 nmpunent 

Material 

- 

•o 

, -13 

ftjilo  eaui 

Thickness 

irm 

Active  Pump 
length  (cm . 

Otame  ter 
i mm  i 

Polarization 

I «lin»r 

BK-? 

1.  507 

134 

l.o 

L 

SF-7  (kitiMJt  leoe 

BK-T 

I 507 

I M 

1 0 

C 

SF-6  Input  lene 

BK-7 

1.507 

1 24 

1.0 

C 

*0  mm  Ampl 

EV-2 

ISO 

0.05 

37 

30 

50 

C 

SF-7  (aitput  I <*n* 

BK-7 

1 507 

1 24 

1.0 

C 

SF-7  Input  I ena 

BK-7 

1. 507 

1. 24 

1.0 

C 

Nirlwr 

BK-7 

1.507 

1 24 

l.o 

L 

W5-  Uirvfci* 

BK-7 

1.507 

1 24 

l.o 

L 

PCS-tryatal 

KD*P 

1.466 

1 00 

10.0 

L 

PCS-  Window 

BK-7 

1.507 

1 24 

l.o 

L 

polarizer 

BK-7 

1.507 

1 24 

l.o 

L 

SK-4  (Jutput  1 eoa 

BK-7 

1.  507 

1 24 

l.o 

C 

SF-6  ln»iut  I en* 

BK-7 

1.507 

1.24 

l.o 

c 

Mill  Ampl 

FV-2 

1.50 

0.  95 

3? 

30 

64 

c 

Beamaplttter 

BK-7 

1.507 

1 24 

l.o 

1. 

Beamsplitter 

BK-7 

1.507 

1 24 

l.o 

L 

SF-5  Output  I-ena 

BK-7 

1.507 

1 24 

l.o 

c 

SF-5  Input  Lena 

BK-7 

1.507 

1.24 

l.o 

c 

6411  Ampl 

g-66 

1.53 

1.22 

37 

30 

64 

c 

Beamsplitter 

*BK-7 

1.507 

1.  24 

1.0 

L 

Heamaplltte  r 

BK-7 

l.  507 

1. 24 

1.0 

L 

SF-4  Output  lens 

BK-7 

1.507 

1 24 

l.o 

L 

SF-4  Input  lens 

BK-7 

1.507 

1 24 

1.0 

L 

Polarl  ter 

BK-7 

1.507 

1 24 

1.0 

L 

Faraday  Rut. 

FR-4 

1 556 

1.55 

2 54 

L 

Polarizer 

BK-7 

1.507 

1.24 

1.0 

L 

641  Ampl 

Q-B6 

1.53 

1.22 

37 

30 

64 

C 

SF-3  Output  lens 

BK-7 

1.507 

1. 24 

1.0 

C 

SF-3  Input  Lena 

BK-7 

1.507 

1. 24 

1.0 

C 

Polarizer 

BK-7 

1.507 

1.24 

1.0 

L 

PC -4  Window 

BK-7 

1.507 

1 24 

1.27 

L 

PC -4  Cryatal 

KD*P 

1.468 

1.00 

6.2 

L 

PC  *4  Window 

BK-7 

1.507 

1 24 

1.27 

L 

Polarizer 

BK-7 

1.50 

1.24 

1.0 

L 

40-mm  Ampl 

EV-2 

1.50 

0.55 

37 

30 

40 

C 

SF-2  Output  Lena 

BK-7 

1.507 

1 24 

1.00 

C 

SF-2  Input  Lena 

BK-7 

1.507 

1. 24 

1 00 

C 

3o-mm  Ampl 

EV-2 

1.50 

0.55 

37 

30 

30 

C 

8F-1  Output  Lena 

BK-7 

1.507 

1.24 

1.0 

C 

SF-1  Input  lens 

BK-7 

1.  507 

1 24 

1.0 

C 

PC3  • POL’* 

- 

16 -mm  Ampl 

EV-2 

1.50 

0.55 

37 

30 

L 
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(Kornd  Energy 
lienalty  Profile 


0.  2S(0.  25  • I)  199  r*i 


0.  25l'i.  36  • 0.173  r2-44, 


0.1240.3*9  • n.  173  r2  sl 


0.22(0  3*5  . 0.  173  r2' 5I 


0.  25(0.  45  • 0.05  r2> 


0.25(0.5  . 0.13*  r2'2,l 


0.  25(0.80  . 0.3  r2  - 0.05  r4l 


TABLE  2 


INPUT  ENERGY  VS  OUTPUT  ENERGY 


Input 

(pJ> 


50  ps 

Output  Fig. 

(J)  Ref. 


1.707 

39.  19 

13 

1. 575 

39.61 

14 

1.313 

38.  28 

15 

1. 182 

36.34 

16 

1.050 

33.69 

17 

0.788 

27.01 

18 

0.525 

19. 13 

19 

0.  263 

10.20 

20 

0. 131 

5.28 

21 

100  ps 

5.25 

73.98 

23 

4.73 

74.75 

24 

4.  20 

74.24 

25 

3. 15 

70.90 

26 

2.36 

61.12 

27 

1.84 

51.87 

28 

1.58 

46.51 

29 

1.31 

40.64 

30 

1.05 

34.  19 

31 

TABLE  2 


INPUT  ENERGY  VS  OUTPUT  ENERGY  (Cont) 


200  ps 


Input  Output  Fig. 

(pJ)  (J)  Ref. 


21.00 

144.77 

34 

18.38 

141.84 

35 

15.75 

139. 12 

36 

13.65 

136.47 

37 

12.60 

134.65 

38 

11.55 

132.28 

39 

10.50 

129.20 

40 

9.45 

125.28 

41 

7.88 

117.47 

42 

5.25 

98.01 

43 

1.31 

40.65 

44 

1 ns 

262  . 56 

272.96 

46 

223. 18 

267  . 77 

47 

170.66 

258.76 

48 

131. 28 

249.43 

49 

105.02 

241.07 

50 

78.77 

229.73 

51 

52.51 

212.66 

52 

26.26 

180.63 

53 

13.  13 

145.62 

54 

2.63 

66.02 

55 

1.31 

40.65 

56 
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